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This  Memorandum  is  a  part  of  RAND's  continuing  VELA  Analysis  study 
for  the  Advanced  Research  Projects  Agency  and  deals  with  a  technique 
for  detecting  nuclear  tests  in  the  atmosphere.  It  presents  the  basic 
physical  models  for  the  infrasonic  effects  which  have  been  observed 
from  nuclear  bursts.  It  discusses  the  propagation  of  low-frequency 
(less  than  one-tenth  of  a  cycle  per  second)  sound  from  nuclear  explosions, 
attempting  to  predict  the  type  of  signal  observed.  For  long-period 
waves,  the  Memorandum  represents  a  refinement  of  previous  work;  for 
short-period  waves,  the  Memorandum  is  perhaps  the  first  serious  examina¬ 
tion  of  the  problem  in  this  context.  This  study  should  be  of  use  in 
revealing  more  effective  methods  of  data  treatment  to  better  evaluate 
the  infrasonic  method  of  detecting  nuclear  explosions. 

The  Memorandum  was  also  submitted  for  publication  in  a  special 
issue  of  the  IEEE  on  high-altitude  nuclear  detection. 

The  author,  a  Consultant  to  The  RAND  Corporation,  is  a  professor 
in  the  Department  of  Aeronautics  and  Engineering  Mechanics,  University 
of  Minnesota. 
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SUMMARY 

This  Memorandum  deals  primarily  with  the  effects  of  geometric 
dispersion  on  low-frequency  mechanical  waves  generated  by  nuclear 
explosions.  This  dispersion  is  the  result  of  the  stratification  of 
the  atmosphere  (to  be  distinguished  from  dispersion  due  to  changes  in 
physical  characteristics  due  to  changing  frequency).  It  is  found 
that  the  pressure  signal  can  be  divided  naturally  into  an  early- 
arriving  acoustic-gravity  portion  (treated  in  this  Memorandum)  and  a 
later--by  about  5  p'r  cent  of  the  travel  time--acoustic  portion. 

In  general,  both  portions  of  the  signal  are  inversely  proportional 

to  the  range  (geometric  spreading  included)  although  at  very  great 

1/6 

ranges  portions  of  the  gravity  wave  fall  off  faster  by  r  ;  the  ef- 
feet  of  dispersion  is  to  reduce  the  signal  by  r  2.  Most  of  the  signal 
is  composed  of  many  propagating  modes  which,  at  a  given  time  and  range, 
will  each  demonstrate  a  characteristic  frequency.  A  simplified  treat¬ 
ment  of  this  complicated  model  picture  is  presented  here.  It  is 
argued  that  a  ray  treatment  for  the  higher-frequency  portion  is 
appropriate.  It  is  shown  that  the  frequency  of  the  long-period  signal 
increases  with  time.  So  long  as  the  frequency  of  the  received  signal 
is  less  than  the  characteristic  frequency  of  the  initiating  explosive 
impulse,  it  is  found  that  the  signal  has  a  universal  form  for  the 
fundamental  mode.  Such  characteristics  as  yield,  range,  and  phase 
velocity  merely  change  the  scaling  of  the  signal.  It  is  concluded 
that  attenuation  is  probably  not  important  for  the  low-frequency 
signals  (below  one  ''ps)  usually  observed. 
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I .  INTRODUCTION 

This  Memorandum  deals  primarily  with  the  effects  of  geometric 
dispersion  on  low-frequency  mechanical  waves  generated  by  nuclear 
explosions.  Geometric  dispersion  in  general  reduces  the  signal 

-k  ( 1) 

amplitude  by  an  amount  proportional  to  r  ,  where  r  is  the  range. 

That  this  is  so  can  be  seen  as  follows:  Suppose  the  amplitude  of 
the  wave  tram  is  A.  The  time  duration  of  the  disturbance  after 
some  time  becomes  proportional  to  the  range  as  a  result  of  dis¬ 
persion.  Total  energy  must  be  conserved  (neglecting  attenuation) , 
so  we  have  approximately 

2 

A  r  =  constant 

Hence  A  ~  r  ,  as  suggested.  This  reduction  is  in  addition  to  that 
caused  by  geometric  spreading. 

The  characteristics  of  low-frequency,  atmospheric,  acoustic- 
gravity  waves  have  been  investigated  by  a  number  of  authors.  Such 
waves  can  be  generated  by  natural  events  as  well  as  by  nuclear  explo¬ 
sions.  We  shall  restrict  the  discussion  largely  to  the  latter.  No 
attempt  will  be  made  here  to  compile  a  complete  bibliography  on  the 
subject.  However,  we  should  mention  normal  mode  calculations  per¬ 
formed  by  Pfeffer,^  Pfeffer  and  Zarichny,^  Harkrider^  and 
Pierce. ^  Harkrider  has  succeeded  in  calculating  the  expected  low- 
frequency  signal  (hereafter  this  term,  and  the  term  long-period, 
mean  waves  with  periods  of  more  than  one  minute) .  He  used  five  modes 
in  the  rather  complex  calculations.  It  is  the  major  purpose  of  this 
Memorandum  to  attempt  to  obtain  the  essential  characteristics  of  the  low- 
frequency  signal  more  simply  (and,  of  course,  less  accurately). 
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The  appropriate  mathematical  theory  for  the  long-period  character¬ 
istics  of  acoustic-gravity  waves  is  given  in  Ref.  2.  From  the  outset 
we  distinguish  between  those  waves  for  which  gravitational  effects  are 
important  (periods  greater  than  60  sec)  and  acoustic  waves  (periods 
less  than  about  60  sec).  For  the  former,  a  normal  mode  treatment 
is  indicated--as  discussed  in  Refs.  2  through  5.  We  examine  some ‘of  the 
expected  characteristics  of  such  waves  on  the  basis  of  simplified 
models.  For  intermediate  periods  (from  30  to  120  sec)  the  situation 
is  one  where  we  must  still  use  normal  mode  methods,  but  where  a 
number  of  normal  modes  can  propagate.  For  the  acoustic  part  of  the 
signal,  a  very  great  number  of  modes  can  propagate  and  a  ray  treat¬ 
ment  is  more  useful. 

We  introduce  here  a  discussion  of  the  normal  mode  treatment. 

For  our  purpose  it  will  be  assumed  that  we  observe  the  mode  in  the 

radiation  zone,  where  kr  »  1,  with  k  a  typical  propagation  constant. 

til 

Then  we  have  ter  the  pressure  of  the  n  mode 

ik  ((o)r-iiut 

P  =  G(r)E  (cu.Z  )  E  (u>,Z)  e  n  (1) 

n  'n  o'n 

where  G(r)  is  a  geometric  spreading  factor  with 

G(r)  =  r"^  (2) 

for  propagation  in  a  flat-earth  atmosphere  (see  discussion  below) , 

■A 

E  is  the  excitation  function  for  point  excitation,  7  is  source 
n  o 

height,  Z  is  receiver  height  (set  equal  to  zero  for  ground-level 

For  details  of  the  excitation  process,  see  Ref.  4. 
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detection),  k  =  cu/c  ,  with  c  the  phase  velocity  of  the  ntn  mode, 
n  pn  pn 

It  is  assumed  in  Eq.  (1)  that  the  excitation  is  a  single-frequency 
continuous  wave.  For  the  usual  pulse  excitation,  we  take  the  Fourier 
Transform  (F.T.)  of  the  initiating  disturbance  to  find 


P 

n 


■  ^  J 


E  (oi,Z  ) 
n  ’  o 


En(«u>Z) 


ik  (co)r-io3t 
n 

e  do> 


(3) 


where  we  assume  the  excitation  pressure  p(t)  has  the  F.T.  P 


P(a>)  =  J*  ela>t  p(t)  dt  ,  (4) 

The  function  p(t)  is  found  by  taking  the  pressure  at  a  range  far 
enough  from  burst  so  that  the  overpressure  is  a  small  fraction  of  the 
ambient  pressure,  but  not  r  far  that  the  inhomogeneous  characteristics 
of  the  medium  influence  the  pulse  form.  Admittedly,  it  may  be  difficult 
to  find  such  an  intermediate  range  for  large  yield  events.  In  such  a 
case,  we  mean  by  p  the  linearized  pressure  disturbance  which  would  be 
observed  in  an  infinite  homogeneous  medium.  The  empirically  observed 
form  of  the  pressure  is  approximately^ 


p(t)  =0  t  <  0 

t  't/,;+ 

P(t)  "  PQ(1  "  ~)  e  t  >  0 


(5) 
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where  t  ®  Yi/-i  t. ,  t,  1.0  sec  and  Y  is  the  yield  in  KT.  The  over- 
+  11 

1/3 

pressure  p  is  also  proportional  to  Y  .It  is  emphasized  that  the 
numerical  value  of  t  is  somewhat  uncertain.  The  F.T.  of  Eq.  (5)  is 

?(o))  =  p  - —  (6) 

-1  2 

(t+  "  ^ 

We  discuss  the  calculations  which  are.  available  for  the  functions 
appearing  in  Eq.  (3).  For  simplicity  the  calculations  have  been  per¬ 
formed  assuming  a  horizontal  atmosphere.  This  is  adequate  for  ranges 
which  are  less  than  a  quarter  of  the  circumference  of  the  earth.  For 
greater  ranges  we  should  take  into  account  the  lack  of  further  geo¬ 
metric  spreading,  due  to  the  spherical  earth.  It  seems  reasonable  to 
include,  approximately,  such  effects  by  modifying  the  relation,  Eq.  (2), 
for  the  geometric  factor.  Thus,  for  the  usual  case  where  wavelengths 
are  small  compared  with  Rg ,  the  radius  of  the  earth,  we  have 


G(r)  =  ]R  sin  (7) 

e 

with  r  measured  along  the  surface  of  the  earth.  The  focus  at  the 
antipode,  r  =  jt  Rg,  is  smeared  due  to  the  fact  that  different  earth 
paths  have  different  temperatures  and  winds  and  thus  different  propa¬ 
gation  times.  The  maximum  value  at  the  antipode  is 


There  are  two,  largely  unresolved,  difficulties  in  the  determination 
of  t+.  First,  it  might  be  that  with  a  large  burst,  the  first  impulse 
takes  on  a  nonspherical  form  before  it  becomes  linearized,  due  to  atmos¬ 
pheric  inhomogeneities.  In  such  a  case  it  is  necessary  to  use  a  source, 
function  more  complicated  than  a  simple  source  function.  Second,  the 
time  of  the  positive  phase,  t  ,  is  not  well  known.  In  Ref.  6,  p.  117,  the 
value  of  t_^_  has  not  attained  its  asymptotic  value  at  the  maximum  range 
shown.  On  the  basis  of  known  information,  this  value  must  be  estimated 
(the  author  is  indebted  to  T.  F.  Burke  of  the  RAND  Corporation  for  this 
statement) . 


where  a ,  order  of  a  few  per  cent,  is  the  fractional  change  in  group 
velocity  for  different  paths.  Other  effects  aside,  we  see  the  focusing 
eff< ct  might  be  appreciable.  An  interesting  consequence  of  Eq.  (7)  is 
that  after  traveling  a  fraction  of  the  circumference  of  the  earth,  the 
signal  shows  no  further  reduction  in  amplitude  due  to  geometric  spread¬ 
ing.  Observed  reduction  in  signal  amplitude  with  range  then  must  be 
due  to  possible  attenuation  (usually  negligible  in  situations  of 
greatest  practical  interest)  or  dispersion,  which  is  the  main  effect. 

In  order  lo  correctly  calculate  the  excitation  functions  and  the 
modal  phase  velocities,  Cpni  we  must  specify  the  atmospheric  charac¬ 
teristics  as  a  function  of  height.  For  this  purpose  Pfeffer  and 
(3) 

Zarichny  used  the  COSFAR  standard  atmosphere  shown  in  Fig.  1.  For  a 
COSPAR  atmosphere  to  300  km  (constant  thereafter)  they  find  the  compli¬ 
cated  phase  and  group  velocity  curves  shown  in  Fig.  2  for  the  first 
few  modes.  In  the  same  figure  the  results  for  a  52-km  atmosphere 
(constant  above  52  km)  are  shown.  It  is  seen  that  this  simpler  model 
gives  all  of  the  characteristics  of  the  300-kra  model  if  attention  is 
restricted  to  the  nearly  horizontal  portions  of  the  dispersion  curves. 
This  is  true  although  the  52-km  fundamental,  for  instance,  is  composed 
of  portions  of  dispersion  curves  from  several  of  the  300-km  modes.  It 
is  seen  that  the  major  portion  of  the  atmosphere  (below  52  km)  deter¬ 
mines  the  propagation  characteristics.  For  periods  greater  than 
60  sec,  it  is  seen  that  the  300-km  atmosphere  group  velocity  curves 
are  either  nearly  horizontal  or  nearly  vertical.  The  near-vertical 


6 


100  200  300  400  500  600  700  800  900  1000  MOO  1200  m  sec  1 
TEMPERATURE  *A 

249  99.5  224  398  S22  895  12)9  1592  20)5  2488  3011  ^583 


Fig.  1— Average  vertical  variation  of  absolute  temperature 
in  the  atmosphere  (estimated  by  the  Comittee  on  Space 
Research;  details  up  to  120  km  given  in  insert.  From  Ref.2) 
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300  km  COSPAR 


Fig. 2— Phase  and  group  velocity  curves  for 
COSPAR  model  atmospheres  (From  Ref. 2) 
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portions  contribute  very  little  to  the  propagating  signal.  Consequently, 
there  is  at  any  given  period  only  one  principal  propagating  mode  (and 
this  one  is  very  similar  to  the  fundamental  mode  for  the  52-km  atmos¬ 
phere).  Of  course,  the  propagating  mode  in  fact  comes  from  different 
modes  of  the  300-km  atmosphere.  For  peritds  from  somewhat  above  30  sec 
to  less  than  60  sec  we  see  that  there  are  other  contributing  modes  in 
the  300-km  atmosphere.  The  important  ones  correspond  to  the  funda¬ 
mental  ana  first  acoustic  modes  for  the  52-km  atmosphere  (there  are 
presumably  more  than  these  two  modes,  although  they  were  not  calculated). 
As  suggested  above,  these  shorter-period  modes  (say,  with  periods  less 
than  10  sec)  are  best  treated  as  rays- -which  in  fact  are  made  up  of  a 
great  many  propagating  modes.  It  should  be  pointed  out  that  there  may 
be  doubts  about  the  300-km  calculation  discussed  here.  Tolstoy^  has 
shown  that  there  may  be  acoustic  modes  trapped  at  about  100  km  between 
the  upper-altitude  temperature  increase  and  the  lower- altitude  density 
increase.  It  is  expected  that  such  an  effect  will  not  influence  the 
propagation  of  periods  of  interest  here  because  of  the  attenuation  at 
such  altitudes. 

It  is  interesting  to  note  that  the  52-km  fundamental  mode  is  made 
up  of  acoustic  modes  from  the  300-km  calculation  for  periods  less  than 
250  sec  and  of  gravity  modes  from  the  300- kra  calculation  for  greater 
periods. 

We  shall  in  what  follows  generally  restrict  the  discussion  to  the 
simpler  dispersion  curves  taken  from  the  calculations  for  thp  52-km 
atmosphere.  We  see  that  tha  long-period  gravitational  contributions 
might  be  expected  to  arrive  wich  a  group  velocity  about  5  per  cent 
greater  (5  per  cent'  less  travel  time)  than  the  shorter-period  acoustic 
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signals.  Of  course,  this  observation  applies  to  the  quasihorizontal 
portions  of  the  dispersion  curves.  The  reader  is  referred  to  Refs.  2, 
3  and  6  for  more  details  concerning  the  normal  mode  calculations. 

In  the  next  section  we  shall  simplify  the  problem  by  considering 
the  behavior  of  a  given  normal  mode  at  great  ranges.  It  will  be  seen 
that,  as  sketched  above,  due  to  the  dispersion  caused  by  the  strati¬ 
fied  inhomogeneity  of  the  atmosphere,  the  signal  amplitude  reduces 
with  range.  This  and  other  characteristics  of  the  signal  will  be 
obtained  from  Eq.  (3).  The  discussion  of  possible  attenuation  effects 
appears  in  Section  IV. 
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II.  SIGNAL  DISPERSION 


For  long-period,  acoustic-gravity  waves,  we  use  the  method  of 
stationary  phase  to  approximately  evaluate  Eq.  (3).  For  large  r 
and  t,  P  and  are  slowly  varying  functions  of  cu  compared  with  the 
exponential.  Then,  following  standard  methods,  we  look  for  the  station¬ 
ary  point  of  the  argument  of  the  exponential,  which  is  seen  to  be  uj 
given  by 


k' 

n 


ton) 


(9) 


Here  (o^  is  the  n  mode  frequency,  if  there  is  one,  whose  group  velo¬ 
city  is  equal  to  r/t.  Then  we  have  for  that  stationary  contribution, 
at  the  frequency  ^--given  by  Eq.  (9) 


ik  (uj  )r  -  iu)  t  _ 

Pn  =  GCr)e  "  "  n  En(VZo)  Sn(VZ)  P(»n>  X 


X  c 


2-rrrc 


2 


gn 


(10) 


using 


k"  = 
n 


'gn  /  =  d 

1  C  £  *  11  Q 

,2  ’  gn  do)  gn 

'gn 


(ID 


obtained  from  Eq.  (9).  Equation  (10)  is  valid  in  general  when  c1 
is  not  too  small  (see  condition  Eq.  (15))  and  when  the  range  is 
large,  usually 


gn 


k  r  »  1 
n 


(12) 


cwm1— ii  ■ 


'  ^ 


T 


li 


The  observed  signal  may  be  a  sum  of  many  modal  contributions  of  the 
type  given  in  Eq.  (10).  If  c^-*0,  as  for  instance  when  the  period 
greatly  exceeds  60  sec,  we  must  modify  the  treatment.  This  is  done 
next . 


We  note  from  Fig.  2  that  for  longer  periods  the  52-km  fundamental 
phase  and  group  velocities  approach  approximately  constant  values. 

It  is  desirable  to  calculate  the  expected  behavior  for  this  situation. 
The  first  acoustic  mode  from  the  52-km  atmosphere  can  be  treated  using 
Eq.  (10).  We  assume  for  long  periods  (a)  «  o)q)  for  the  fundamental 


mode ,  n  =  0 


k>) 

o 


Ac  3 

_ R 

2  2 

c  Oi 
peo  o 


(13) 


with  (u  the  transition  frequency  (about  2tt/t  where  T  w  60  sec)  and 
o  o  o 


-  Ac  (-^) 
P 


poo 


(14) 


where  c  is  the  phase  velocity  for  very  long  periods  and  Ac  is  of 
p00  P 

the  order  of  the  change  in  the  phase  velocity,  to  be  obtained  by  fitting 


the  phase  velocity  of  the  fundamental  52-km  mode  in  Fig.  2.  We  see 

Ac 

from  that  figure  that  — 0.05. 


The  restriction  on  c'  for  the  validity  of  Eq.  (10)  can  be  written 

gn 

conveniently  in  terms  of  Eq.  (14).  It  is 


\3/ 2  /  .  x  1/2 
6  \  u?  /  \  r  ]  \  Ac. 


1/2 


«  1 


(13) 


It  might  be  thought  from  a  casual  observation  of  Fig.  2  that  there 
is  a  similar  situation  for  short  periods.  However,  when  the  curves 
are  put  on  a  linear  rather  than  a  logarithmic  pJot,  it  is  found  that 
they  actually  have  finite  slope  as  the  period  tends  to  zero. 
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with  X  the  wavelength  corresponding  to  c  .  For  values  taken 

o00  '  *  ?- 

from  Fig.  2,  the  left  clue  of  Eq.  (15)  is  .02  for  a  range  of  8000  km; 
for  most  ranges  of  interest  Eq.  (15)  is  fulfilled. 

From  Eqs.  (5)  and  (6)  note  that  for  ordinary  yields  Y,  the  con¬ 
dition  that  to  «  to  implies  that  we  need  the  low-frequency  limit 
behavior  of  P  in  Eq.  (3).  From  Eq.  (6)  we  see 

P0t2(-io)>,  tt)  «  t'1  (16) 

Similarly,  it  is  reasonable  to  believe  (and  calculations  in  Ref.  4 

show)  that  the  characteristic  frequencies  of  E  are  of  order  to  . 

o  o 

Then  we  also  need  the  long-period  behavior  of  E^,  seen  to  be  approxi- 

k 

mately  constant;  this  value  of  E  is  called  E  below.  Substituting 

o  Lr 

these  low-frequency  approximations  in  Eq .  (3)  we  find 


Po,LP  ‘  G«  %'+  ELP(Z) 


ELP^Zo^  dt 


(17) 


with 


*LP  3  t  J  exp 


Ac  3 

no  ,  _ R 

c  2  2 

p00  c  to 
.  p®  O  j 


r  -  <ot  ^  dtu 

J 


(18) 


*LP  resPonse  to  3  -sita  function  excitation)  is  found  by  stand¬ 


ard  methods  to  be 


k 

See  p.  5309  of  Ref.  4. 
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1  ~1_  ,  r 

3  Vi  ,t  ''c 


-1/2 

c)|  V2!/ 

3 


>>r- 

poo 


-  t  >  0 


in  r  ->*'2  (  r  ,  3/2 

TTst (t  ■  ^->]  {J1 <2  [f  <*=  -  r~>]  > 

pco  J  w  ~  U  poo  -J 


V + j  i a  \* (t  -  f->j/2>}  ^  -  c-  > 

\  ~ -x  P=°  J  p=° 


with 


Ac  1/3 

T  =  (“jj2,  *^>  (20) 

c  cu 

poo  O 

Here  an(^  ^1/3  are  t^le  uSual  Bessel  and  modified  Bessel  functions. 

It  is  easily  seen  from  Eq.  (18)  that  T  ^  is  the  order  of  the  value  of 
u)  for  which 


Ac  3 
_ E  0L. 

2  2 
C  U3 
p®  O 


that  is,  where  the  cubic  term  in  the  exponential  begins  to  bring  about 
convergence  of  th^  integral.  The  error  made  by  using  the  low- frequency 
approximation  for  P(<u)  in  the  region  to  >  t^1  is  negligible  if 


Ac  riu 
P  ° 


«  1 


This  condition  is  satisfied  for  ranges  exceeding  a  hundred  kilometers 
The  result  of  Eqs.  (17)  through  (20)  shows  that  the  medium  is  essenti 
ally  a  low-pass  filter  for  long-period  waves  when  there  is  concern 
only  with  the  most  important  mode,  the  fundamental. 


15 


III.  DISCUSSION  OF  LONG-  AND  INTERMEDIATE-PERIOD  SIGNAL  DISPERSION 


The  reader  is  reminded  that  here  long  periods  are  those  over  60 
sec  and  intermediate  periods  are  those  between  about  30  and  120  sec. 
Under  the  assumption  that  for  long  periods  the  fundamental  mode  has 
group  and  phase  velocities  which  are  nearly  constant,  the  result 
given  in  Eqs .  (17)- (20)  was  obtained.  The  following  characteristics 
appear  from  Eqs.  (17) -(20)  for  the  early  long-period  component--some- 
times  called  the  gravity  wave: 

1.  The  signal  arrival  time  is  approximately  fixed  by  the  long- 
period  phase  velocity,  c  .  The  precursor  at  earlier  times 
increases  in  amplitude  exponentially.  The  signal  amplitude 
then,  after  this  "crescendo,"  increases  like  J  ,  where 
At  =  t  -  ^ — •  This  increase  in  amplitude  continues  until 

pCO 

>  3-2 

At  ~  3t  t  ,  when  frequencies  are  such  that  the  approxima¬ 
tion,  Eq.  (16),  fails.  For  later  times  the  more  general 
result  given  in  Eq.  (10)  could  be  used  (see  below).  Since 
the  growth  in  signal  amplitude  is  due  to  the  increase  of  P(u>) 
with  (ju,  an  alternative  procedure  can  be  used.  We  use  the 


pulse  transform,  Eq.  (6),  and  the  period  as  a  function  of 
—  given  in  Eq.  (22)  to  approximately  correct  for  the  pulse 

T 

form.  For  this  purpose,  multiply  Eq.  (19)  by 


t.  ,2  .  -1 


M-a  a 


(21a) 


We  assume  in  such  a  treatment  that  the  approximate  effect  of 
the  pulse  spectrum  is  to  multiply  the  signal  at  time  At  by 
the  value  of  che  spectrum  function  appropriate  to  the  frequency 
occurring  at  that  time. 


f 


f 
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2.  The  first  important  signal  arrival,  when  t  -*  r/c  ,  ^as  a 

long  period  followed  by  successively  shorter-period  arrivals 

(normal  dispersion).  From  Eq.  (15)  it  is  seen  that  using  Ttn 

for  the  approximate  zeros  of  the  Bessel  functions,  when 
3/2 

(At/r)  >  1,  the  period  for  signal  delayed  by  At  is 


We  see  the  period  decra.tsnr;  like  the  ^  power  of  the  delay  in 
this  normally  dispersed  part  of  the  signal. 

3.  The  range  dependence  is  interesting.  Tha  amplitude  of  the 
early  part  of  the  low-frequency  signal  is  seen  from  Eqs .  (17)  ~ 
(20)  to  fall  off  like  the  2/3  power  of  the  range  due  to  dis¬ 
persion  (remembering  to  take  the  derivative  as  in  Eq.  (17)). 
For  ranges  considerably  less  than  ,  we  use  Fq.  (2)  for  the 

function  G  and  find  that  the  long-period  wave  shculd  decrease 

7*  /  6  ^ 

in  amplitude  like  r  .  For  greater  ranges,  r  ~  R  ,  the 

-2/7 

dependence  is  roughly  like  r  “.  The  period  or  the  first 

% 

arrival,  at  a  given  At,  increases  like  r  according  to 
Eq.  (22).  As  the  range  increases,  the  amplitude  decreases 
and  the  period  increases  due  to  dispersion.  Relatively 
slow  decreases  in  amplitude  with  range  have  been  found 
experimentally  for  signals  which  circle  the  globe. ^  ^ 

4,  From  Eq.  (17)  and  the  scaling  laws  discussed  in  Section  I# 
we  see  that  the  amplitude  is  proportional  to  yield.  The 
signal  waveform  is  otherwise  independent  of  yield  so  long 
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3-2 

as  At  <  3t  t  .  This  should  be  true  for  the  long-period 
+ 

gravity  wave,  but  care  must  be  used  when  the  delay  is  such 
that  the  wave  merges  with  the  rhorter-period  components. 

3-2  -1 

5.  For  At  <  3t  t  ,  so  that  tu  <  t,  ,  we  see  the  signal  has  a 

+  -r 

universal  form.  The  yield,  range,  and  atmospheric  effects 
enter  the  scale  of  the  amplitude  (for  yield)  and  th.e  time 
scale  (for  range  and  atmospheric  effects) . 

These  characteristics  are  also  obtainable  using  the  simple  station¬ 
ary  phase  treatment  given  by  Eq.  (10).  It  is  emphasized  that  this  re¬ 
sult  applies  only  to  the  fundamental  mode  (as  defined  for  the  52-km 
atmosphere).  Often  one  or  more  of  the  higher  modes  can  be  expected  to 

appear  also,  and  this  effect  must  be  calculated  using  Sqs,  (9)-(10). 

4  4c 

Consider  an  example:  r  '•*  lO  km,  u)q  *  2tt/60,  =  .05,  c^  ^ 

poo 

-km/sec.  Then  we  find  t  «  50  sec.  Accordingly,  the  result,  Eqs. 

3 

(17)- (20) ,  is  valid  for  At  <  3r3  t+2;  and  with  Y  =  4  MT,  t+  -  16  sec, 

3 

we  need  At  <  24  min.  From  Eq.  (22)  for  At  *  10  sec,  the  period, 

T,  is  115  sec;  for  At  *  4000,  T  ^  55  sec. 

The  universal  result,  Eqs.  (17)-(20),  is  shown  in  Fis.  3.  Thera 
the  quantity  t2  “  ILp  (~~)  is  plotted.  The  characteristics,  dis¬ 
cussed  above,  can  be  seen.  The  dashed  envelope  curve  gives  the 
approximate  modification  due  tc  a  4*fIT  pulse  form,  given  by  the  factor 
of  Eq.  (21a)  (parameters  for  this  curve  are  those  of  the  example  just 
completed).  It  is  interesting  to  note  that  the  initiation  of  the 
signal  using  the  velocity  c  is  correctly  messured  from  the  first 
positive  maximum.  Ths  precursor  (At  <  C)  is  not  actually  in  conflict 


F»g. 3 — Low -frequency  universal  response  of  the  atmosphere  for  the  fundament,  mode 


t2  I  .«(At/r) 


•  i  t 

O  -K  NJ  o  ro  o 


8T 


I 


19 

with  the  general  theorems concerning  the  velocity  of  propagation  of 
the  signal  wava  frant  bec&'ise  of  the  form  of  the  dispersion  curve.  We 
should  use  Eqs.  (S}-(10)  to  calculate  the  curve  when  ~  is  3uch  that 
0)  2:  u)^;  for  the  example  above  (4  MX)  this  means  for  periods  of  about 
one  in  the  dimensionless  units  of  the  graph.  However,  the  approximation, 
Eq.  (13),  is  not  bad  for  much  higher  frequencies,  so  that  it  is  not  un¬ 
reasonable  to  accept  the  results  for  greater  time  delay. 

We  turn  now  to  a  discussion  of  the  intermediate  periods  (periods 
between  30  and  120  sec).  We  use  Eq.  (10)  and  in  general,  we  must  con¬ 
sider  many  propagating  modes  in  the  signal.  The  number  of  modes  may  be 
quite  large  when  the  period  is  much  less  than  60  sec.  We  see  the  fol¬ 
lowing  characteristics  for  the  intermediate  component: 

* 

1.  By  referring  tc  Ref.  4,  it  is  found  that  the  fundamental  mode 
for  periods  greater  than  90  sec  has  an  excitation  function 
which  is  constant  at  an  altitude  of  18,5  km.  More  work  needs 
to  be  done  on  the  excitation  functions  in  order  that  the  fre¬ 
quency  dependence  may  be  examined  more  clonclv.  As  already 
remarked,  we  assume  that  for  periods  greater  than  60  sec  the 
excitation  functions  are  constant. 

2.  Other  things  being  equal,  the  amplitude  of  the  n*”*1  mode  is 
inversely  proportional  to  the  square  root  of  the  slope  of  the 
modal  group  velocity  for  the  frequency  in  question.  As  sug¬ 
gested  above,  the  most  important  modes  are  those  whose  group 
velocity  curves  are  most  nearly  horizontal  (see  Eq„  (10)). 

3.  The  general  characteristics  of  the  signal  are  determined  as 
follows  (see  sketch  below  wnich  is  a  schematic  representa- 

See  p.  5309  of  Ref.  4. 
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ticn  ol  the  first  three  modes).  First,  for  a  given  range 
and  time,  Eq.  (?)  xs  used  to  determine  the  appropriate  group 
velocity.  Then  the  lower  cuives  for  52-km  atmosphere  are 
used  from  Fig.  2  to  find  which  periods  appear  for  the  various 
modes.  Unfortunately,  the  calculations  given  there  are  only 
for  the  fundamental  and  the  first  acoustic  modes.  The  other 
acoustic  modes  must  be  imagined  as  lying  above  the  first. 

The  slope  of  the  group  velocity  curve  is  determined  from 
the  calculation.  The  result  is  that  for  any  range  and  time 
there  can  be  a  great  number  of  modes  arriving,  each  with  its 
appropriate  period  as  explained.  We  see  in  the  sketch  the 
three  periods  for  the  first  three  acoustic  modes  at  r,t. 
Others  are  not  shown.  At  a  given  time  many  different  modes, 
each  with  its  own  frequency,  will  make  up  the  signal,  though 
the  fundamental  mode  is  dominant  for  low  frequencies. 


••CCJ* 
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4.  The  periods  of  the  various  modes  decrease  as  the  time  of 

observation  increases  (normal  dispersion) ,  as  can  be  seen  from 
the  sketch.  Normal  dispersion  is  (usually)  observed  experi¬ 
mentally.  The  characteristic  empirical  period, 

2rro)o\  is  shown  in  the  sketch.  The  period  for  which  the 
pulse  form  becomes  important,  2rrt+,  is  also  shown;  it  is  the 
period  below  which  uhe  simple  form  of  the  fundamental,  Eqs . 
(17)-(20),  must  be  modified  because  of  pulse  form. 
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The  rue st  prominent  frequency  in  the  intermediate  signal  cannot 

be  determined  without  knowledge  of  the  excitation  functions. 

For  short  periods  where  a  ray  treatment  applies,  we  may  expect 

the  signal  to  show  a  frequency  maximum  at  that  of  the  pulse 

maximum,  u)  ^  **  t+*.  Hence,  the  "frequency"  of  the  signal 

1/3 

would  be  proportional  tc  Y  '  . 


7.  From  Eqs .  (6)  and  (10)  it  is  expected  that  the  amplitude  of 

the  signal  would  be  proportional  to  pQ,  that  is,  proportional 
1/3 

to  Y  .  Howe.er,  the  excitation  function  can  indirectly, 
through  the  parameter  t+  in  P,  influence  the  dependence  of  the 
amplitude  upon  yield.  The  acoustic  signal  amplitude  can,  of 
course,  be  influenced  by  fluctuating  atmospheric  propagation 
characteristics . 

Before  closing  this  section,  it  is  useful  to  express  the  inter¬ 
mediate-period  mode  results  in  terms  of  the  power  spectrum  function. 


We  use  the  definition  of  the  power  spectrum  for  a  segment,  2Tq,  of 


record  f(t) 


P .  S .  (u>)  =  jp-  |  Aj  (uj) 


with 


^  (uj)  =  •—  J  °  e  iUJt  f(t)  dt. 


The  result  for  the  power  spectrum  of  the  n  ’  mode  is  found  from 
Eqs.  (10)  and  (24)  to  be 


*»***»?*  J 
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where  K  is  constant,  independent  of  range,  and  9  is  the  azimuth  angle 
about  the  source.  The  result  must  be  independent  of  the  range  for  a 
loss- free  medium.  The  total  energy  is  given  by  a  sum  over  terms  like 
Eq.  (27),  covering  all  of  the  propagating  modes;  we  neglect  inter¬ 
ference  effects  which  do  not  affect  the  energy  content.  The  range 
independence  is,  of  course,  a  reflection  of  the  fact  that  geometric 
dispersion  in  the  medium  does  not  affect  the  energy  of  the  propagating 
wave.  This  suggests  using  the  integral  of  the  square  of  the  observed 

pressure  signal  since  this  quantity  would  show  no  reduction  with  range 

2 

due  to  dispersion.  It  would  decrease  like  1/r  rather  than  1/r  .  The 
range  dependence  of  the  energy  would  be  entirely  geometric.  For 
r  >  Re,  the  average  geometric  reduction  should  be  slight.  For  such 


ranges,  we  should  find 


constant,  independent  of  r 


(28) 


f  p2(t)  dt  - 
0 

Noise  can  be  removed  by  subtracting  an  integral  like  Eq .  (28)  for  the 

noise  signal.  It  should  be  emphasized  that  Eq.  (28),  representing  as 

it  does  a  fraction  of  the  explosion  energy,  should  be  proportional  to 

yield.  This  dependence  should  be  compared  with  the  fractional  power 

yield  dependence  of  some  other  characteristics.  We  see  from  Eq.  (28), 

as  remarked  above,  that  since  dispersion  causes  the  signal  to  spread 

in  time  proportional  to  range,  we  can  expect  the  pressure  to  be  propor- 

tional  to  r  for  ranges  such  that  r  >  . 

For  shorter-period  acoustic  waves,  we  can  employ  the  full 

formalism  of  ray  treatment.  This  will  not  be  done  here.  However, 

a  few  outstanding  characteristics  are  worth  noting.  First,  the 

range  dependence  due  to  atmospheric  stratification  is  also  r  2 

as  for  normal  modes.  We  might  expect  interference  effects  due  to 

multi-path  phenomena  for  shorter  ranges.  The  most  prominent  circular 

+  -1 

frequency  in  this  part  of  the  signal  will  occur  at  (t  )  ,  the 

characteristic  frequency  of  the  initiating  explosion. 


:»■ 
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IV.  SIGNAL  ATTENUATION 

We  shall  consider  in  this  section  a  few  of  the  more  prominent 
attenuation  mechanisms.  F?rst,  we  examine  absorption  due  to  molecular, 
vibrational  relaxation.  Using  accepted  vaLues  we  have  for  the 
attenuation,  coefficient 

q  2 

a  ~  10'y  — — 2~2~  cm"  ,  T»  .002  -  .05  sec  (29) 

1  +  (1)  t 

2tt  _ 

For  0)  =  a  5-sec  period 

C*  ~  10  11  cm  1  (30) 

(L 

The  signal  reduces  to  a  fraction  of  its  value  in  10 '  km  due  to  molecular 
absorption.  For  periods  above  a  few  seconds  we  conclude  that  such 
attenuation  can  be  neglected. 

Consider  now  the  attenuation  which  would  be  caused  by  scattering 
from  velocity  irregularities  set  up  by  atmospheric  turbulence.  Due 
to  the  pulsed  nature  of  the  initiating  signal,  we  shall  for  simplicity 
assume  that  once  a  signal  is  scattered  it  is  lost.  It  is  true  that 
some  of  such  energy  is  regained,  but  for  the  present  purpose  of  estima¬ 
tion  we  neglect  such  contributions.  Since  we  are  interested  here  in 
long-period  disturbances  whose  wavelengths  are  usually  greater  than 
the  scale  of  atmospheric  irregularities,  we  restrict  the  discussion 
to  Rayleigh  scattering.  We  find  for  the  attenuation  coefficient 

a  ~  (ka)4  -  (AN) 2 ,  ka  <  1  (31) 

(AN)2,  ka  ~  1  (32) 

(Eq.  (32)  is  given  for  completeness.) 
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In  Eq.  (31),  a  is  the  turbulence  scale  iengch--the  correlation  length, 

AN  is  the  turbulent  velocity  fluctuation  relative  to  the  speed  of 
sound,  and  we  have  assumed  that  the  turbulent  eddies  are  "close- 
packed,"  that  is,  the  medium  is  fully  irregular.  The  Rayleigh  result 
goes  over  to  a  geometrical  optics  result  of  the  order  given  in  Eq.  (32) 
for  higher  frequencies.  A  typical  scale  length  for  atmospheric 
turbulence  is  about  1  km.  Taking  the  maxi.miur  attenuation  given  in 
Eq.  (32),  for  a  period  less  than  about  20  sec,  and  using  AN  =»  ,005, 
a  value  appropriate  to  the  upper  atmosphere,  we  find 

O'  *  ~  40,000  km 

Longer  periods  would  give  even  less  attenuation.  We  conclude  that 
because  of  the  low  contrast  of  turbulent  velocity  fluctuations, 
attenuation  due  to  such  scattering  is  negligible  for  low-frequency  waves. 

For  higher  frequencies,  although  the  scattering  from  turbulent 
eddies  is  slight,  there  can  be  appreciable  reflection  from  laminar 
wind  ducts. 

More  generally,  it  can  be  argued  from  experimental  observations 
(see,  for  instance,  Refs.  8  and  9)  that  there  is  no  appreciable  signal 
attenuation  at  present  ranges  of  observation.  Any  attenuation  mech^  'sm 
would  give  a  certain  reduction  in  signal  per  unit  length.  Then,  if  the 
Signal  is  reduced  by  attenuation  to  a  fraction  by  traveling  a  fraction 

of  the  way  around  the  earth,  it  would  be  reduced  by  that  fraction  to 

th 

the  n  power  cir.a  path  n  times  as  long.  But  the  first  reduction 
due  to  attenuation  must  be  negligible  since  no  such  drastic  reduction 
is  observed  in  subsequent  passes  around  the  earth. 


r. 


V.  CONCLUSIONS 


We  have  examined  the  effect  of  dispersion  on  acoustic-gravity 
waves  generated  oy  nuclear  events  at  great  ranges.  It  is  shown  by 
using  the  normal  mode  calculation''  of  Pfeffer  that  it  is  natural  to 
divide  the  signal  into  three  parts:  low-frequency  gravity  modes, 
intermediate-frequency  signals,  and  acoustic  signals.  The  charac¬ 
teristics  of  the  first  two  were  examined  in  Sections  II  and  III.  It 
is  found  that,  in  general,  dispersion  causes  a  reduction  of  signal 
amplitude  proportional  to  the  square  root  of  the  range,  regardless 
of  the  portion  of  the  signal  considered.  This  is  in  addition  to  the 
reduction  due  to  simple  geometric  spreading.  The  additional  range 
dependence  is  sufficient  to  give  the  observed  reduction  in  signal 
level  over  that  expect.ed  from  ordinary  cylindrical  spreading. 

We  summarize  the  results.  In  order  to  calculate  the  pressure 
in  different  period  ranges  (alternatively  for  different  time  delays), 
see  below. 


Fundamental  Mode 


For  a  period  greater  than  —  =  60  sec  and  2P  t+  (alternatively 
Ac  he 

At  less  than  3  — r  and  3  — ^  use  Eqs.  (17)-(20), 

c  c  L  w  Z  tZ 

pro  p®  O  + 


where  it  is  assumed  the  excitation  function  is  constant  in 


frequency. 

For  a  period  -less  than  either  above  (alternatively  At  greater 
than  corresponding  quantities  above),  use  Eqs.  (9)-(10). 


o 


For  all  period  use  Eqs.  (9)-(10). 


o  o 


It  is  appropriate  to  comment  on  the  expected  variability  of 
pressure  records.  First,  it  seems  reasonable  to  assume  that  the  gravity 
waves  will  exhibit  little  variability  with  changes  in  tropospheric 
conditions,  since  they  have  very  great  wavelengths.  The  same  cannot 
be  said  of  the  acoustic  waves  (periods  less  than  about  20  sec)  .  Their 
wavelengths  are  of  the  order  of,  or  less  than,  the  size  of  prominent 
features  of  the  lower  atmosphere.  From  Fig.  1  it  is  seen  that  there 
is  an  acoustic  velocity  maximum  at  about  50  km.  This  maximum  has  a 
width  of  about  10  km,  which  corresponds  tc  the  wavelength  cf  a  30-sec 
period  wave .  For  such  periods  and  shorter,  we  should  expect  that  changes 
in  wind  velocity  at  50  km  might  seriously  affect  the  modal  excitation 
functions  and  to  a  lesser  extent  the  dispersion  curves,  or  alternatively 
the  lay  propagation  characteristics.  In  this  way,  the  large  and  variable 
winds  at  such  altitudes  might  cause  considerable  variability  in  the 
acoustic  and  perhaps  intermediate-frequency  signals.  This  would  be 
particularly  true  since  relatively  small  changes  in  wind  speed  are 
seen  from  Fig.  1  to  be  sufficient  to  cause  the  disappearance  of  the 
lower  duct  altogether,  as  observed  from  the  ground. 

Another  conclusion  which  should  be  emphasized  is  that  the  integral 
of  the  square  of  pressure  record,  the  "energy"  in  the  wave ,  may  serve 
as  a  more  sensitive  characteristic  for  detection  than  the  signal 
itself.  This  is  so  since  such  a  measure  removes  the  effect  of  disper¬ 
sion.  Thus,  tne  range  dependence  is  just  that  of  the  (range  periodic) 
geometric  spreading.  It  is  seen  that  this  energy  is  proportional  ro 


vie  Id . 


For  future  wci'k  it  is  suggested  that  power  spectral  functions 
(P.S.F.'s)  be  calculated  for  the  sigral  records  in  order  to  deter¬ 
mine  the  signal  characteristics  along  the  lines  discussed  in  this 
Memorandum, 

Signal  variability  forms  an  interesting  topic  in  itself.  It 
would  be  helpful  to  attempt  to  correlate  this  variability  with  upper 
atmosphere  winds.  It  will  be  interesting  to  examine  further  P.S.F.'s 
to  see  the  behavior  of  thr  spectral  peaks.  In  particular,  the  reduction 
in  frequency  with  increased  time  of  travel  should  be  looked  for.  Such 
frequency-travel  time  results  can  bs  used  to  obtain  an  experimental 
determination  of  the  group  velocity  dispersion  curves  for  the  various 
modes.  Finally,  it  seems  desirable  to  obtain  the  characteristics  of 
the  modal  excitation  functions  for  shorter  periods  and  as  a  function 
of  altitude.  It  is  noted  that  these  functions  may  yield  important 
information  on  burst  height  (see  Eq.  (1)). 

The  results  of  the  .implified  calculations  discussed  in  this 

(4) 

report  have  been  compared  with  the  more  complete  results  of  Harkrider, 
The  frequency  dispersion  and  signal  amplitudes  correspond  well,  with¬ 
in  10-20  per  cent  in  the  cases  considered. 
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